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Abstract 

Following the general principles of noncommutative geometry, it is possible to define a metric 
on the space of pure states of the noncommutative algebra generated by the coordinates. This 
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give some examples of computations in both commutative and noncommutative cases. 
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I Introduction. 



Though particle experiments are going further in energy and consequently deeper in the 
structure of matter, the geometric structure of space time is still unknown. Classical differential 
geometry does not allow to take seriously into account both general relativity and quantum 
mechanics since the latest renounces intuitive geometric concepts while the first grounds its 
decription of gravitation on purely geometric concepts. Noncommutative differential geometry! 
gives a mathematical framework for a geometric understanding of fundamental interactions. 
Saying geometric understanding, one would like to say clearer understanding: for instance the 
noncommutative standard model!'!! gives a geometric interpretation of the Higgs field together 
with an estimation on the mass of the corresponding boson. 

But it is still difficult to draw an intuitive picture of a noncommutative space, as one can 
do for an euclidean or even riemannian space. A noncommutative space is described by a C*- 
algebra A, a faithful representation of A over a Hilbert space TC, and an operator D acting 
on TC. D has a compact resolvant and is possibly unbounded. To be precise, the algebra is 
restricted to the norm closure of the set of elements a G A such that [D, vr(a)] is bounded. 

A distance is then defined on S(A), state space of A, by 

d($,tf) = sup{|$(a)-tf(a)| / || [A a] II < 1}, V$,f G5(i). (1) 

a<=A 

If A is commutative, pures states correspond to characters. Thanks to Gelfand construction, 
they are interpreted as points, and A is the algebra of functions over these points. When A is 
not commutative, this interpetation is no more possible but the distance formula between pure 
states remains unchanged. 

There is not always a clear physical interpretation of this distance. When A is the algebra of 
smooth functions over a riemannian spin manifold, TC the space of L 2 -spinors and D the classical 
Dirac operator, then the noncommutative distance coincides with the geodesic distance. When 
A is tensorised by an internal algebra, for instance the diagonal 2x2 matrices, then one obtains 
a space of two sheets with geodesic distance over each of them. Noncommutative distances 
have been consideredHi'i in the case of finite algebras. Moreover, the classical distance in one 
dimensional lattices can be obtained via the noncommutative approach!'!. 

In Refs. |iO|JTT| , the problem is introduced in a general framework: let L be a Lipschitz 
seminorm on a partially ordered real vector space A. L determines a metric pi in the state 
space S(A): 

Pl{p, v) = sup { |/i(a) — v(a)\ / L(a) < 1} with p,, v 6 S(A). 
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Pl determines a Lipschitz seminorm L PL over the space Af[S(A)] of affine functions on S(A): 
L PL (f)= snp { li^zjpi /v^v} with f e Af[S(A)). 

As A is isomorphic to a dense subspace of Af[S(A)], the question is: under which conditions 
one has L PL = L ? Answers are given in Ref. [11]. In the noncommutative framework, taking 
L(a) = || [D, a] ||, the question becomes: how to characterize the metrics p coming from a Dirac 
operator ? 

In this paper, we carry out the calculation of distances in spaces associated to finite di- 
mensional algebras. It seems natural to restrict to finite dimensional representation of these 
algebras since, in the noncommutative approach of the standard model, the internal space is the 
space of fermions (more mathematical arguments to avoid infinite representations of finite di- 
mensional algebras can be found in Ref. |l|). Therefore, A is a direct sum of k matrix algebras, 
since any involutive algebra over C which admits a faithful finite dimensional representation in 
a Hilbert space is a direct sum of matrix algebras. For k — 1, the simpliest interesting case is 
A = M2(C). The associated space is a fiber space whose base has only one point and the fiber 
is C 2 . For k = 2, we study the noncommutative space associated to A = M P (C) © C, p G N. 
This is a two point space with fiber C p over one of the point. Some applications can be found 



in Ref. [13], where A = M 2 (C) © C is used to build a first model of quantized spacetime. 

For k > 3, we restrict to commutative algebras. Then A = C k and S(A) is simply a set of 
k points. We choose 7i = <C k . For the three point space with any real selfadjoint operator D 
and the n point space with some particular operators D, we explicitly compute distances. To 
find a Dirac which gives a desired metric, it is enough to inverse formula. This is not possible 
in the four point case for we show that generic distances are roots of polynomials which cannot 
be solved by radicals. A possible solution consists in modifying our definition of commutative 
spaces. Using a slightly more complicated representation of A over a space 7i' larger than 
Ti,, one shows that there always exists an operator D' giving the desired distances. Moreover 
(A, T~C', D') is a real spectral triple which fulfills all the axioms of noncommutative geometry. 



II General result and notations. 

All along this paper, A is a unital C*-algebra represented in a Hilbert space Ti. D is a 
selfadjoint operator on 7i which does not belong to the commutant of A. Its components are 
Dij, 1 < i,j < n. S(A) is the set of pure states of A. The distance associated to the triplet 
(A, H, D) by ([1|) is denoted by d. 
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Lemma 1. Let A + be the subset of positive elements of A, and $, \& G S(A). 



d($,^) = sup { |$(a) - tf(a)| / ||[D,o]|| = 1}. 

ae.4 + 



Proof. Let = arg (($ — \l/)(ao)) where Oq E A reaches the supremum in ([]]), namely: 

\\[D,ao]\\ < 1 and |($ - *)(a )| = dist($, *). 

The supremum is also reached for the selfadjoint element bo = ^(a e~ ie + a^e 19 ) G A since 

|| [A&o] ||<I^ + MU <!, 



The same is true for Co = &o + H^oll-^ ^ *4+> so we restrict to A + . 
Suppose now ||[D,co]|| < 1. Take e = ||p° ]|| £ then 

||[Ae ]|| = l and |$(e ) - *(e Q )| = > l*(o>) - *(<*)!, 

which is impossible since Co is chosen to reach the supremum. So ||[D, co]|| = 1. 

If the supremum is not reached, the proof uses a sequence {a n } of elements of AM 

• Once for all, any element a £ A that appears in a proof is selfadjoint. 

• The canonical basis of C n ( or M. n in case) is denoted by |1), |2), \ n). 

• When A = M n (C), a pure state is determined by a normalized vector £ G C n : o;^(a) = 
£*a£, Va G M n (C). Two normalized vectors determine the same pure state if and only if they 
are equal up to a phase. In other terms, S(M n (C)) = CP n_1 . 

• For any unitary operator U of H, the jauge transformed of is u^(a) = ui^iUaU^ 1 ), 
Va G A. If D = U~ 1 DU, then df)(u>£, cJ^) = d D {uj^ u<*) for 

sup{|(^-^ c )(a)|/||f/- 1 J Df/,a]|| = l} = sup {|(^ — ^)(f/af/ _1 )| / || [L>, f/af/ -1 ] || = 1} . 
ae.A UaU- x eA 



III One point space. 

The first non trivial example with a single matrix algebra is A = M 2 (C), represented in 
H = C 2 bjA3a = ( au ai2 V D is a 2 x 2 selfadjoint matrix with two real and distinct 



a 21 a 22 

eigenvalues D\,D2. S(A) = CP 1 is isomorphic to the sphere S 2 . An explicit one to one 
correspondence is the Hopf fibrationB: £ = (£1,^2) £ CP 1 is associated to (a^,6^,c^) G S 2 by 

a 5 = 2TZe{^ 2 ), 6 5 = 21m(^ 2 ) and c e = l^ 2 - |6| 2 - 



Proposition 2. 



2yi-i(e,oi 2 

«H^c) = \D 1 -D 2 \ — ^ C « = C C' 
= +00 otherwise. 



Proof. Let C/ be the unitary operator of 7i such that D = U 1 DU = diag(-Di, D 2 ). A direct 
computation yields ||[D, a]|| = |«i2 1 l-Ci ~ D 2 \, thus the norm condition in ([[]) becomes 

2 

Furthermore, (u^ — u^)(a) = £* a £ — £*a£ = £ ay(&£j — GO)- 

j,j = l 

^ I 7^ |Ci|) then a with all coefficients zero, except an = L, verifies the norm condition 
(|) for any L G K + , and |(c^ -w c )(a)| = Z-||6| 2 - |Ci| 2 |. Thus d 5 (u; 6 a; c ) = +00. 
If |6| = ICil, then |6| = K2I since ||f || = ||C|| = 1. So c € = c c and 

\(u 6 - u c )(a)\ = \27Ze (a 12 (66 - C1C2)) I < 2|a 12 ||(&6 - CiC 2 )|. (3) 

As any vector of CP™ -1 is defined up to a phase, we assume that £1 = Ci is real. Let 9% = arg(£ 2 ) 
and 9^ = arg(C 2 ). Take an = a 22 = and arg(a 2 i) = — % — 9^). Then 

\Ke (a 12 (e"i6 - C1C2)) I = |a 12 | yCi |6l |2sin(^=-^)| = |a 12 ||66 - CiCal- 

a reaches the upper bound in and verifies the norm condition (0) as far as one chooses 
I «i2 1 = | Di ^ g2 | ■ So df)(uj£,UJ(;) = \ Dl l D2 \ [£i£ 2 — C1C2 1 ■ An intrinsic formulation is found by 
writing 

tt* _ cc* _ [ _ £i£ 2 — C1C2 
^ ^ ~ V 66 - C2C1 

Tr(a* - CC*) 2 = 21^6 - C1C2I 2 . Developing the trace yields 1 - |<£, 0| 2 = |66 - OGIY thus 

dfi(^,wc) = i^-i^ v^-lfcOl 2 - 

Finally, df)(uj^,u)^) = d^iU" 1 ^, U~ X C,) which is equal to (^(c^, c<^) as seen before. ■ 

We say that two states oj^uj^ G S(A) are at the same altitude if = c^. By an easy 
calculation, for two such states, d(u^LU^) = ^ t D2 \ y/ (Qg ~ <\) 2 + ~ ^c) 2 - m other terms, 
up to a constant factor, d is nothing but the euclidean distance restricted to plans of constant 
altitude. The distance between to planes of different altitudes is infinite. 

In a one point space with a fiber of higher dimension than C 2 , one needs an explicit formula 
for the norm of a selfadjoint n x n complex matrix. These cases are not studied here. 
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IV Two point space. 

Consider the algebra A = M n (C) © C represented on H = C n © C by A 3 a 



x 
y 



( m \ 

a; G M n (C), y E C The simpliest interesting operator is .D = ( „ q J> with m G C n a 
non zero column vector. S(A) is the union of the single pure state of C, u = Identity, with 
S(M n (C)). 

Proposition 3. 

d(u>£,UQ) = u^jj- i/ £ and m are colinear, 
d(uj^,LUo) = +oo otherwise. 

d(uj£,u>(;) = +oo otherwise. 

Proof. We may assume that ||m|| = 1, for dividing D by ||m|| means multiplying distances by 
||m||. Thus, there is a unitary operator u G M n (C) such that m = u\l). With U = ^ J 
l) = U~ 1 DU and z = x — yl n (selfadjoint by lemma 1) one has 

_ n 

ll[Aa]|| 2 = El^r, (4) 
i=i 



since 

— 2r|l> \ ||2 _ || / z\l)(l\z* 
(l\z J 11 ~ 11 I (l\zz*\l) 



\\[D,a] 
and 



n n 

\\z\l)(l\z*\\ = \\z\l)\\ 2 = E \~ Zil \\ \(l\zz*\l)\ = S |^| 2 . 

2=1 1 = 1 

Note that || [D, a] || 2 = 1 implies \za\ < 1 VzG{l,n}. 

n _ 

^ (w { - w ) (a) = C< - 2/ = C*^ = . s 

Assume 7^ for k G {2, n}. Take the matrix z with all coefficients zero except Zkk — L G 
By (f|), 2; satisfies the norm condition of formula ([I]) and |(o>£ — Uq)(o)\ = \ik\ 2 L. Thus 
df)(uj^uj^) = +00. 

Assume £j = 0, Vz G {2,n}: there is constant 6* such that £ = e ld \l). So |(c<j£ — co> )(a)| = 
I -2n I < 1- This upper bound is reached by z with all coefficients zero except Zn — 1. 

n _ _ 

(w{ - w f )(a) = . £ - CO)- 

Assume — CpCz) 7^ f° r P> ^ ^ {2,n}. The proof is similar as i) with z = except 
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Assume — Q(j = Vi,j G {2,n}. This is equivalent to £j = Qe 10 , with # a constant. In 
other words, (£ — (e ld ) ~ |1). Furthermore, since ||£|| 2 = ||C|| 2 = 1, one has |£i| = \d\. Thus 

n _ _ n _ _ n _ _ 

IV ~ ^0)1 = 2 I E ^ e M&i - CiCO) I < 2| - C<Ci)| 

i,j=l 1=2 1=2 



< 2JE\za\\XUi - OCil 2 < 2 J S - OCil 2 - (5) 

n _ _ x 

Take z = except = — 0Ci|( ^ — C*Ci| 2 ) 5 - z reaches the upper bound (0) and 

i=i 

nrz : 

verifies the norm conditions. Thus djj((jJe,uA = 2W E — GCi| 2 - As in proposition 2, an in- 

y i=2 

trinsic expression is found by developing Tr(££* — CC*) 2 ? an d the formulas for dp are the same ■ 



The cases A = M P (C) ©M g (C) is not studied here, neither is the space associated to a sum of 
three or more algebras with at least a noncommutative one. We focus on sums of commutative 

k 

algebras. Then A is isomorphic to ® C. The space associated to k — 1 has no interest. With 

i=l 

k = 2, there is only one distance to compute which is equal to r^r ■ The generic cases k = 3, 4 
with a real operator D, and some examples with k = n G N are considered below. Before, we 
present general results on commutative finite spaces. 



V Commutative finite spaces. 



A n point commutative space is determined by a triplet (A,T~t,D) in which A = ®yC is 

i 

represent ated over Ti = C n as a diagonal matrix: 



A3 a 



( ai ... \ 

: a 2 : 



V o 



where G C, but for distance computing we restrict to a.; G 1R + thanks to lemma 1. To 
make computations easier, we only consider real operators D. As D only appears through its 
commutator [D, a], we assume that it has the following form: 



D 



( D 12 ... 

i ^23 

V #m 







A»-l,» / 



with G 
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Pure states can be interpreted as points of a n point space whose function algebra is A: a(i) = a%. 
The distance between two points i,j of this finite space is 



d(i,j) = sup { |a$ 



I \\[D,a]\\=l}. 



(6) 



In finite spaces, D may be interpreted as the adjacency matrix of a latticei: two points i 
and j are connected if and only if ^ 0. For instance, in the four point space, the restriction 
to -D13 = D 2 4 = corresponds to a cyclic graph: 



A path jij is a sequence of p distinct points (i, i 2 , i p -i, j) with D ikik+1 ^ 0, Vfc G {l,p—l}. 
Since d(l, 2) = in a two point space, the length of 7^ is by definition 

p-i 1 



fe=1 lA^ fc+1 | 



Two points i,j are said connected if there exists at least one path 7^. The geodesic distance 
Lij is by definition the length of the shortest path connecting i and j. 

Proposition 4. 

i) Let D' be the operator obtained by cancelling one or several lines and the corresponding 
columns of D. Then d^ > dr>. 

ii) The distance between two points i and j depends only on the matrix elements correspond- 
ing to points situated on a path joining i and j . 

iii) The distance between two points is finite if and only if they are connected. 

Proof. 

i) Let define e £ A by = if the i th line and column are cancelled, e% = 1 elsewhere, e is 
a projection, commuting with A, and D' = eDe. Thus ||[D',a]|| < ||[D,a]||,Va £ A. So 

sup{|aj — a,j\ j \\[D, a] || < 1} > sup{|aj — aj\ j || [D', a]|| < 1}. 

ii) Let Tij denote the graph associated to the set of points belonging to any path 7^, and 
lij the set of points which are not on any path 7^. Any point of 1^ is connected by at most 
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one path to T^. In other terms, V/ G hj there is at most one point mi G T^- such that I and m; 
are connected and 7; m; has all its points (except mi) in J^. 

Let D' be the operator obtained by cancellation of all lines and columns associated to points 
of Iij. Assume that the supremum in (|6|) with D' is reached by a' G A = C n . Consider a G C™ 
such that a p = aj, except for the points of 1^ for which a; = a m; or a\ = if m/ does not exist. 
Then | \[D, a] \ \ = \ \[D', a'}\\, so do{hj) > dD'(i,j)- By i), d^ii,]) < d' D (i,j). Hence the result. 

Hi) Suppose % and j are connected. There is at least one path 7^ = (1,12, ■ ■-,i P -i,j) whose 

length is the geodesic distance L^. Let obtain D' by cancelling all lines and columns which do 

not correspond to points of 7^-. Then d(i,j) < dD'(i,j)- By the triangular inequality, one has 

p-i p-i 1 

d D '{i,j) < £ d D >(i k ,i k+1 ) = £ — 1 = Lij. 

k=1 k = l \Vi k i k+1 \ 

d(i,j) is smaller than the geodesic distance, thus it is finite. 

If i and j are not connected, define a G C n by Oj = i > 0, a k = a% if and z are connected, 
0^ = otherwise. Then [D, a] = and |a$ — | = t. As t is arbitrary, d(i,j) is infinite. ■ 



• For simplification purpose, we write 

Ojj = — di and 2 = 021 , x« = 2 < i < n — 1. (7) 

In the n = 3 and n = 4 point case, this reduces to 

y = 031 and z = a 4 i, (8) 

• In the four point space, the components of D are 

di — ~^r~ ' ^ 2 = ~FT~ ' = 7T - ' ^ 4 = IT - ' ^ 5 = ~FT~ ' ^ 6 = 7T - • (^) 

-^12 -^13 -^14 -^23 -^24 -^34 



VI Distance on a regular space. 

A n point commutative space is called regular when all coefficients of operator D are equal: 

D = {Dij} = {k(l - Sij)} , keR. 

Proposition 5. 

i) The distance between two points i,j of a regular space of constant k is 
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ii) If the link, and only this link, between two point i\,i 2 is cut, D ili2 = 0, then 

^ 2) = mvS- 



Proof. In the regular space, the problem is symmetrical: all distances are equal and we 
compute d(l, 2). When a link is cut, we take i\ — 1, — 2 to fix notations, and denote by D' 
the operator D with Dyi = . In both case, (H) and (|7|) yield 



d(l,2) = sup{ |x|/||[D or D',a]\\ = 1}. 

aeA+ 

We first compute the norm of the commutator, and then find the supremum. 



(10) 



Lemma 6. 

i) II [A II 2 



n n 

|fc| 2 £ . £ a?, = \k\ 



i=l j=i+l 





n-1 / 


\k\ 2 


x 2 + £ 

i=2 \ 



n-1 



j=i+l 



I „ll|2 



\[D',a 



\_kl_ 

~Y 

\k\ 2 

c 

+ 



n n 



£ . >: ,aj - J(E .>: ^a?.) 2 - 4a? 2 £ o £ a 



"(1,2) 



i=l j'=i+l 
-(1,2) 

n-1 



i=3 j=i+l 



n—1 t . } 

£ I x ~\~ ( x x t ) ~\~ £ ( x^ x j 

i=2 \ i=«+l 



n—1 / n—1 \ n—1 n—1 

( £ ( x? + (x - X;) 2 + £ (xi - x ? -) 2 I ) 2 — 4x 2 E E ( Xi - x 7 -) 2 

1=2 V j=i+l I i=2 j=i+l 



ii) For the regular space, the supremum of x in (1C) is reached when all Xj 's are equal. 



Proof, i) C = i[D, a] is the n x n matrix 

/ ia 12 

C = k 



ia 2 i ' • taij 



\ 



/ 



with rank < 2 since its kernel is generated by the (n — 2) independant vectors 

Afc = (— — ; — 0; 1; 0; ), 1 beeing at the k th position, 3 < k < n. 

0-21 0-21 

Moreover C is hermitian and traceless, so it has two non zero real eigenvalues ±A. Thus 
A = \/ Tr ^ - ■ A direct computation yields 



n n 

A = k\ I E E (an) 2 , 
v i= ij=i+i K %3! 
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n / n—1 n— 1 n— 1 

Finally || [£>,«] II = = \M = |&| A /.E _S (o^) 2 = % 2 + E xf + E _E (a* - x 3 -) 2 - 



Let C = i[D', a]. C" is the n x n matrix 



C = k 



( za 13 \ 
via. 



za 3 i 



V 



/ 



with rank < 4 since ker(C) is generated by the (n — 4) independant vectors 

Ap = (0; 0; 0; 1; 0), 1 beeing at the p i/l position, 5 < p < n. 

CI43 a.43 

Beeing hermitian and traceless, C has four real eigenvalues ±A 1; ±A' 2 . Thus its characteristic 
polynomial is 



X(C) = x n ~\x 2 - a; 2 )(x 2 - a 2 z ) = x n - (a; z + \' 2 z )x n -' 2 + a; z a 2 z x 



/ 2\ v n-2 i \/2\/2-tA n _4 



IF 



A direct computation yields A' x 2 + A' 2 2 = \Tr(C 



n n 

fc 2 E E (a,,-) 2 . 
i=i j=i+i 

"(1,2) 



The coefficient of X n ~ 4 is the sum of all the minors of C of degree 4. A minor M(l,k,l,p) 
composed with the first (or second column) and three others columns k,l,p ^ {1,2} (and the 
associated lines) is also a minor of C . As C is of rank < 2, its minors of degree greater than 2 
are null, so M(l, k, l,p) = M(2, k, l,p) = 0. The same is true for the minors M(q, k, l,p) with 
q ^ {1, 2}. Finally, the only non zero minors are 



M(l,2,/,p) = k*Det 



( iau iai p \ 

ia 2 i ia 2p 

ian iai2 ia tp 

\ ia pl ia p2 a pi ) 



k'Det ( au aip 

Oj21 «2p 



k CL21O.pl- 



,/ 2\/2 



Summing all these minors gives X[ A 2 = a 12 ^ ^ yu, z 

i=3 p=Z+l 



E E ?/a 2 ; . Then, solving (0) yields 



| [17, a] 



) E E a 2 -+ /(E E a?.) 2 -4a 2 2 E E - 2 

2 \ i=l j=i+l %J i=lj=i+l lJ/ 7=3 j=i+l 

(1,2) V -C 1 - 2 ) 



n—1 n—1 n—1 

Let f(x,X2,---,x n -i) = x 2 + E x 2 + (x — x^) 2 + E E (xj — x,) 2 and suppose that 

i=2 i=2j=i+l 

(x,X2, ...,x n -i) G R n_1 reaches the supremum, namely 



f(x,x 2 , --,x„_i) 



and d(l, 2) = \x\ 



then 
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a) x is positive: from the global parity of /, x can be chosen positive. 

b) Xi < §, Vi G {2, ...,n — 1} : suppose that p of the x^s are greater than | and denote 
them generically by x p . Consider now the (n-l)-uplet in which all Xp S £1X6 replaced by |. Then 
/ decreases for x 2 + (x — x p ) 2 > ^ + (x — |) 2 and (xj — x p ) 2 > (xj — §) 2 . Fixing the values of 
the remaining Xj's leads to see / as a function of the single variable x: 



/(x) = x 2 + p(f f + Ex 2 + p(x - §) 2 + E(x - Xif + Ep(§ - Xi) 2 + EE(x, - x,) 2 , 



» 3 



f(x) = 2x + 2px + 2E(x - Xi) + Ep(§ - x;). 



As Xj < ^ < x , /'(x) > when x > 0. / is continue and lim f(x) = +oo, so there is x > x 

x— »oo 

such that /(xo) = pjp. In other terms, the initial (n-l)-uplet in x p does not reach the supremum 
which is in contradiction with our hypothesis. So p = 0. 

c) Xi > 0, Vz G {2, n — 1}: the proof is the same by replacing all Xj < by |. 

y4Z/ Xi's are equal: let A and A be the two smallest value of the Xj's, with A < A . If 
A = A then it comes immediatly that all X-i S £1X6 equal. If A ^ A then 

< A < A < Xi < | , Vz G {2, ...,n- 1}. 

Assume that m of the Xj's are equal to A. Summing over Xi ^ A, one obtains: 

f(x, x 2 , x n -i) = x 2 + mA 2 + Ex 2 + E(x — x«) 2 + m(x — A) 2 + Em(A — x«) 2 + E(xj — x,-) 2 . 

i i % i,j 

Fix the values of x, 7^ A and consider now A not like a constant but like the value of a variable 
x min . Then / can be seen as a function f m of the two variables x min and x with 

df m 



linn 



As ^^-(x min ,x) < for x min G [A,A[, one has / m (A,x) < / m (A,x) = p^p. Moreover 

-^(A, x) = 2x + 2m(x - A) + E2(x - x;) > 0, 
ax i 

so there is x > x such that f m (x ,A) = -^p, which is inconsistent with our hypothesis. So 

A = AM 

Proof of proposition 5. i) According to lemma 5, Xj = X2, 2 < z < n — 1. The norm 
condition in ( |I0|) becomes 

2( n - 2)x 2 + [2(2 - n)x]x 2 + [(n - l)x 2 - _] = 0, 

l«r 

which has no real solutions in x 2 unless |x| < my~- This upper bound is reached when 



°° 2 2 ~ 2\k\ y n- 
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n— 1 n— 1 n— 1 

m) Let /ii(x,Xj) = £x 2 + (x — Xj) 2 , h (xi) = S S (xj — x,-) 2 , g(x, Xj) = /ii(x, a:,) — 2x 2 . 

i=2 i=2 j=i+l 

Lemma 6 yields 



II [D\ a] f = ffi (/n + /i + yjh\ + hl + 2g.h^j • (12) 
Let xo = sup {x//ii(x, Xj) = j^p}. As (7 and ho are both positive, ( |T2"D implies that 2) < xo- 



Imitating i), one finds this upper bound is reached when all equal, and x = 




In finite spaces which are not regular, distances are not always explicitly computable. The 
cases n = 3 and n = 4 are considered below. 



VII Three point space. 



Proposition 7. For a three point space with operator 



D 



( D 12 D 13 \ 

D 12 ,023 

V D Vi D 23 J 



Dij el, d(l, 2) 



£?3 + £>2 2 3 



F)2 n2 , n2 n2 , n2 r>2 
-^12-^13 1 -^12-^23 1 -^23-^13 

i/ie others distances come from suitable permutations of indices. 



Proof. Equation (f|) and notations (BJ) gives 



(2(1,2) = sup {x / || [D, a] 



-0>i 2 x ~D 13 y 
D 12 x D 23 (x-y) 

Di 3 y D 23 (y-x) 



!}• 



By direct calculation, \\[D, a]\\ = a/ 'D 23 (x — y) 2 + D r3 y 2 + D\ 2 x 2 . Thus (2(1,2) is the largest 
value of x for which there is a point (x, y) belonging to the ellipse 



(D 2 3 + D 2 2 )x 2 + {D\ 3 + D 2 23 )y 2 - 2D 2 23 xy = 1. 



(13) 



(2(1, 2) is the positive x for which the equation in y (13) has a zero discriminant, that is for 



x 



Dfs + ^3 



D 2 2 D 2 13 + D 2 12 D 2 23 + D 2 23 D 2 3 
The three distances verify an inequality of the triangle "power two" since 



rf(l,2) 2 + d(2,3) 2 > d{l,3f 



(14) 
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and two others inequality by permutations. The question rises of inverting the problem: for 
three positive numbers (a,b,c) verifying (II), is there any operator D giving (a,b,c) as non- 
commutative distances ? 



Proposition 8. Let a,b,c G R + verifying a 2 + b 2 > c 2 , b 2 + c 2 > a 2 , a 2 + c 2 > b 2 . Then, 
there is an operator D such that d(l, 2) = a, d(l, 3) = b, d(2, 3) = c, explicitly given by 



D 
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2{b 2 + c 2 -a 2 



{a + b + c)(-a + b + c)(a-b + c){a + b-c) 
Di 3 and -D 23 coming from permutations of a,b,c. 

Proof. Writing jk- = i? i2 , 7^- = -R23, -rt = R13, proposition 7 gives 

U 12 U 23 U 13 

1 1 1 

+ 



rf(l,2) 2 R 12 R23 + R13' 
d(l, 2) 2 is the electrical resistance between points 1 and 2 of the triangle circuit. Finding Di 




means finding three elements Rk p that induce a resistance d(i, j) 2 between points A classical 
result@ indicates that the triangle circuit is equivalent to a stellar circuit (r 1? r 2 ,r 3 ) with 

R12 = — {nr 2 + nr 3 + r 2 r 3 ), (15) 
^3 

i?i3 and -R23 coming from suitable permutations of indices. In the stellar circuit 

d(l, 2) 2 = n + r 2 ( 2n = rf(l,2) 2 + d(l,3) 2 -d(2,3) 

d(l,3) 2 = n +r 3 which implies < 2r 2 = d(l, 2) 2 + d(2, 3) 2 - 3) 
d(2, 3) 2 = r 2 + r 3 [ 2r 3 = d(l, 3) 2 + d(2, 3) 2 - d(l, 2 

Inserting in (|T5| ) this gives 



2(d(l,3) 2 + d(2,3) 2 -d(l,2) 



12 



2(d(l, 2) 2 d(l, 3) 2 + d(l, 2) 2 d(2, 3) 2 + d(l, 3) 2 d(2, 3) 2 ) - d(2, 3) 4 - d(l, 3) 4 - d(l, 2) 4 
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VIII Four point space. 



In a three point space, distances can be computed for any real operator D. This is not true 
for a four point space. Notations and @ are used. 



Theorem 9. 

i) d(l, 2) is a root of a polynomial of degree 5 < 12. 

ii) d(l, 2) is generically not computable: it is a root of a polynomial not solvable by radicals. 

iii) It is computable in the following case: when i = i = oo, 

if d\ < d\, else 
if d\d§ = d^d^, else 
»/ C < 0, else 



d(l,2) 



di 



ri 2 (rf 2 +rf 2 )(rf 2 +rf 2 ) 
((^(it-dide) 2 



max 



d 2 (d|+d 2 ) 

i3+d 4 ) 2 + (di-d 6 ) 2 ' 



'i3-rf4) 2 + (rfl+d6) 2 ^ 



w/iere C = ((d 3 + d 4 ) 2 d 6 + (di - d 6 )(d 3 d 4 - d 6 2 ))((d 3 - d 4 ) 2 d 6 + (di + d 6 ) (d 3 d 4 + d 6 2 )). 



d(l,3) 



a/ d 3 2 + 4 
a/ d 4 2 + d 4 



max 



\J (di d 3 +d 4 d 6 ) 2 



\J (d 1 d 3 +d 4 d 6 ) 2 



lf(d 2 3 + dl)<(d 1 d 6 -d 3 d i ) 2 , 

if (d\ + d\) < (did 6 - d 3 d 4 ) 2 , 
otherwise. 



A/(d 3 +d4) 2 + (rfl-rf 6 ) 2 ' A/(d 3 -d4) 2 + (dl+rf 6 ) 2 ' 

Permuations of the d\s give d(2, 3), d(3, 4), d(l, 4) (Vesp. d(2,4)) from d(l,2) (resp. d(l,3)). 



The rest of this section is devoted to the proof of theorem 9. 



[D,a] 



For "r* 



/ 





X 

di 


y 

d 2 


z 
d 3 


\ 




X 







x—z 






dl 


d 4 


d 5 






y 







y-z 






d-2 


d 4 


de 




\ 


z 
d 3 


z—x 
d 5 


z-y 

d 6 





/ 









(!) 







Va e Aa 



(x, y, z) e M 3 , define the functions: 

„2 „,2 „2 „,\2 



a r 



^_ jr_ + <r_ + (x-y) + ( x - z ) + (y - *0 

' d 3 d 4 d^ dfi 



di d 2 



'6 



g(y - g) gCg - y) vi z - g) 

d^dg d 3 d 4 d 2 ds 

a( r*) + A/a(T*) 2 - 4/5(Y > ) 2 , 
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and the surfaces M and T: 

M = {r ) GK 3 /jj(T ) )=2}, 

T = {7*eR 3 / /( r>) = 0}, with J\f c 

Lemma 10. 

i) ForaeA + , \\[D,a)\\ 2 = 

ii) For G TV suc/i £/ia£ a(~r > ) = 2, grad(/)("r > ) = 0. 



Proof, i) The four eigenvalues of a] are Aj = ±^y a ± \fct 2 — A(3 2 {j2), so 

||[^,a]|| 2 = Vc, 2 -4/3 2 )(r^). 

We show that f^(~^ > ) = 0, the proof beeing the same for the other components of grad(/). 
As T> G J\f C T and a(T > ) = 2, /?(T*) = ±1. If Q{jf) = 1 then 

a(T)=20{l>) and ^L(r>) = — (T>) - 2/3(T>)^(T>) = — (r>) - 2^(T>). 

Explicit calculation of ot("r > ) — 2/3(~r > ) = shows that 

x |/ — z ?/ 2 — x z x — y 



, , , . -( 

9y v / dy 



which leads to ^(T*) — 2§^(~r > ) = 0. The proof is the same if assuming /3(~r* 



Thanks to notations and lemma above, © leads to 

c/(l,2) = sup{(l|? ; :|l)/^GAr}. (16) 
This formula is not useful, for M is not defined by a quadratic form. It is easier to work with T . 

Proposition 11. d(l , 2) G {(1\T>\1) / T* G J 7 and g(T > ) = fj(^) = 0}. 

Proof. The supremum in (|T6| ) is reached at a point "r* such that grad(n)(~r > ), if it is defined, 
is parallel to the x axis. If a(~r > ) = 2, then grad(n)(T^) is not defined but — %(~^) = 

by lemma 10. If a(~r > ) ^ 2, then gra3(/)(T i> ) is colinear to gra3(ra)(T*), so §^(~^) = §^("^0 = 
0. To complete the proof, one just remarks that V~r* G M 3 , there is a G .4.+ such that "r* = r^, 
for instance a = (£, £ — x, £ — y, £ — z) where £ = sup{|x|, |z|}.B 

According to this proposition, the distance is a common root of a polynomial in several 
variables and its various derivatives. Before undertaking explicit calculations, we recall general 
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results about polynomial systems. 



Notes on systems of polynomial equations. 

Let P and Q be two polynomials of the form 

P(x) = a n x n + a n _ x x n ~ l + ••• + a 
Q(x) = b m x m + b m „ l x m ~ 1 + ... + b . 

with a n , b m ^ 0. Without calculating the roots pi, qj of P, Q, one finds by algebraic manipulations^ 
of the coefficients a, and bj the resultant of P and Q: 

Res(P, Q) = a™b n m U( Pi - Qj ), 1 < i < n, 1 < j < m. (17) 

Res(P, Q) is a polynomial in the a^'s and bj's. P and Q have a common root if and only if their 
resultant is zero. A particular resultant is the discriminant: 

Dis(P) = Res(P,P'). 

P has a double root if and only if Dis(P) = 0. If P and Q are polynomials in x, y, z, then 
Res[P, Q, y] denotes the resultant of P and Q seen as polynomials in y. Equivalently Dis[P, z] 
stands for discriminant of P seen as a polynomial in z. 

Proposition 12. Let P(x, y, z) be a polynomial of degree 2 in z, whose coefficients are real 
functions of x and y. If P(x , y , z ) = ^{x , y , z ) = fj(x , Vo, ^o) = for some (x , y , z ) E 
M 3 , then xq is a root of the polynomial Dis[Dis(P, z),y], and yo is a double root of the polynomial 
Dis(P,z)(x ,y). 

Proof. Writing P(x, y, z) = a(x, y)z 2 + b(x, y)z + c(x, y), a direct computation yields 

V = Dis(P,z) = a{Aac-b 2 ), 

dV da , , 9n , db d 9 dc 
— = —{8ac-b 2 )-2ab— + 4a 2 — , 
ay ay ay ay 

„ ,dP OP . , 2 da n , db A ode 
dy dz dy dy dy 

P(x ,y ,z ) = ^(x ,y ,z ) implies V(x ,y ) = 0, i.e. (8ac - b 2 )(x , y ) = b 2 (x ,y ), thus 

dP dP dV 
Res (-g^^^)( x o,yo) = -g-(x ,yo). 
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Therefore, ^(x , y , z ) = q^{x , y , z ) implies -^(x , y ) = = V(x , y ), thus y Q is a double 
root of V(xo, y) and 

Dis{V, y)(x ) = Dis[Dis[P, z],y]{x ) = 0. ■ 



Proof of theorem 9. Proposition 11 and 12 yields 

d(l,2) e W Dis[V{x,y),y] = 0} with V(x,y) = Dts[f{x,y,z),z}. 

Instead of V(x, y), one uses the effective form without the correctif term al^tf^ appearing in 
([[7]), so that zeros of V e ff correspond exactly to the existence of a common roots of / and 

rl Jn 

with /„ the leading coefficient of / seen as a polynomial in z and n = deg(/). Note that the 
numerator is taken after a possible (but not always possible) simplification of the fraction. 

%) By direct computation, V e ff(x, y) = Viy 1 , < i < 4. Exact expressions of the t^'s are 
given in appendix. They are polynomial in x of the form: 

V A (x) = v Ao , V 3 (x) = v 3l x, V 2 (x) = v 22 x 2 + v 2o , Vi(x) = v l3 x 3 + v Xl x, V (x) = Vq^X 4 + v 02 x 2 + v 0o . 



The discriminant J of a polynomial C = Ciy 1 of degree four is 

J(C) = Res[C, C'\ = C^Cs^C^C* 2 - 4d 3 C 3 + lSOAC-A - C (4C 2 3 + 27C C 3 2 )) 

+2(-2C 2 3 (C 1 2 - AC C 2 ) + dC^d 2 - A0C C 2 )C 3 - 3C (^i 2 - 2AC C 2 )C 3 2 )C A 
-(27d 4 - 144C C 1 2 C 2 + 128C 2 C 2 2 + 192C 2 CiC 3 )C 4 2 + 256C 3 C 4 3 ). 

Replacing Cj by Vi(x) shows that J is an even polynomial in x of degree 6 < 12. 

ii) We compute an explicit counterexample: assume d\ = d 2 = d 3 = d A = = 1 and ^- = 0. 
Then 

f(x, y, z) = x 2 + y 2 + z 2 + (x - yf + (y - z) 2 - (x(y - z) + z(x - y)) 2 . 
It is a polynomial of degree 2 in the variables x, y and z. Direct computation gives 

V eff (x,y) = 2-Qy 2 + 3y 4 + 4x 2 (~l + y 2 )-4xy (-l + y 2 ), 
Dis(V eff ,y) = -768 (-54 - 54 x 2 + 135 x 4 + 296 x 6 - 368 x 8 + 128 x 10 ) . 
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Let p(x) = — 128x 5 + 368x 4 — 296x 3 — 135x 2 + 54x + 54. p has one real root x±, two distinct 
complex ones x 2 & X4 and their conjugates x 3 & x$. Galois theory showes that p is not solvable 



by radical (for a comprehensive presentation see Ref. jlTf ): 

p is irreductible over Z because it is irreductible over Z5. Indeed modulo 5 it becomes 
q(x) = 2(x 5 + 4x 4 + 2x 3 + 2x + 2) which has no roots in Z5. Therefore p is irreductible over Q. 

Let E/Q be a splitting field extension of p. As p has five distinct roots, its Galois group 
G = Gal(E/Q) is isomorphic to a subgroup of the symmetry group S5 which is the permutation 
group of X = {x\, X2, x 3 , 24, 25}. As p has no repeated roots, p is separable so |G| = [E/Q] 
where |G| denotes the order of G and [-E7/Q] is the index, i.e. the number of cosets Q in G. If 
a is a root of p then [Q(a), Q] = 5 so |G| = [E/Q] = [E/Q(a)][Q(a), Q] is divisible by 5. Thus 
G contains an element of order 5: the 5-cycle r = (12345). 

The restriction to X of the complex conjugation gives rise to an element o of G: a = 
(23)(45). As a is of order two, |G| is divisible by 2. Moreover to = (124) G G is of order three 
which divides |G|. Thus |G| is a multiple of 5 x 2 x 3 = 30 and divides IS5I = 120. Since 5*5 
has no subgroup of order 30, |G| G {60, 120}. If |G| = 60 then G = A 5 but r £ A 5 . So G = S 5 . 

S n is solvable for n < 4 but is not solvable for n > 5, so G is not solvable. Then, by Galois 
theorem, p is not solvable by radicals. 

in) When d 2 = d 5 = and d^ 6 7^ d^d^. 

Dis{V eff ) = -16^ 6 4 4 4 2 4 4 (^ + dl){x 2 - rf 2 )(x 2 (rf 3 ^4 - <Me) 2 - d\(d\ + dl)(d\ + dl)) 
(x\(d 3 - d A f + K + dsf) - d\{d, - d 4 ) 2 ) 2 (x 2 ((d 3 + d,f + (d, - d 6 ) 2 ) - dj(d 3 + d A ff. 

This polynomial has four single roots ±xo, ±Xi and four double roots ±£2, ±^3: 

a; - |Gti|, xi - \ax\ |d 3 d 4 -rfid 6 | ' 



X2 I "1 1 V ( d3 +d 4 )2+(di-d6) 2 > X3 - l lly (ds-dif + idi+de) 2 - 

By proposition 11 and 12, d(l,2) is one of these Xj's, and the associated yi is a double root 
of V e ff(xo,y) . The corresponding Zi is determined by solving f(xi,yi,z) = 0. Then one 
checks under which conditions each Xi verifies n(xi, yi, z/) = 2 and finally take the greatest one. 
Considering xq, one has yo = di, z = 0, and 



. 1 d\ Hd 2 -diy J 2 if d(<di 

,,(0:0, 2/o,^o) = 1 + ^ + ^ ^ = I 2 | >2 if d?>4 

Therefore xo may be solution only if d\ < cig. Likewise, with the corresponding y\ and zi given 
in appendix, one checks that x\ cannot be solution unless C < 0. On the contrary, x 2 and x 3 
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may always be solutions for there are y 2 , z 2 and 2/3, z 3 such that n(x2, y 2 , z 2 ) = n(xs, y%, z 3 ) = 2 
under no particular condition. By proposition 2, cancelling all links except d\, d(l,2) < xq. 
So if d\ < d\ then d(l,2) = xq. As x\ > x 2 and x\ > X3, d(l,2) = X\ if C < 0, otherwise 
c?(l, 2) = max(xi, x 2 ). 

When dide = d^4, x± is not defined but the proof follws the same way. 

Calculation of d(l, 3) is the same, except we are searching the maximum of y. Dis(V,x) is 
a polynomial in y of degree twelve, with single roots ±yo, ±yi and double roots ±y 2 , ±2/3: 

y = V d 3 2 + d 6 2 , 2/1 = a/ (ii 2 + (i 4 2 , 

Idi rf 3 + c? 4 d 6 | \d1d3 + d^del 

V2 = di = , y 3 = - ■ 

/ 2 2 / 2 2 

y (4 + 4) + (di - 4) y (4 - 4) + (rfi + d 6 ) 

With the associated Xi,Zi given in appendix, one checks that yo (resp. yi) may be solution if 
(g?3 2 + de 2 ) 2 < {d^d^ — did$) 2 (resp. (d\ 2 + c? 4 2 ) 2 < (d 3 d^ — did e ) 2 ). As above, y 2 and 2/3 may 
always be solution. Then, remark that 2/2, 2/3 < 2/o and 2/2,2/3 < 2/1 • Finally, 2/0 and 2/1 cannot be 
simultaneous distinct solutions, for adding both conditions yields 2/0 = 

The four point space shows that there is no hope to find a general formula for the metric in 
any commutative finite spaces: distances cannot be read directly in the Dirac operator through 
a finite algorithm. Computing the metric requires a more pragmatic approach and shall be 
undertaken case by case. 

IX Distances and axioms of noncommutative geometry. 

In the previous discussion, we worked with triplets (A, 7i, D) as if they satisfied all the 
axioms of noncommutative geometry. These axioms are introduced in order to recover the 
standard spin and riemannian geometries in the commutative easel. Accordingly, for our dis- 
tances to be bona fide noncommutative generalizations of riemaniann metrics, they have to be 
computed using triples satisfying all these axioms. 

However, these axioms lead, in the finite cas el, to matrices whose size increases rapidly 
with n and thus prevents any computation except in few simple cases. This is the reason why 
we did not use these axioms up to now, but we shall see that the axioms do not put any con- 
straints on the distances. 

Proposition 13. Let (dij) 1Ki j <ni ^j be any finite sequence of possibly infinite strictly posi- 
tive numbers such that dij = dji and dij < dik + dkj- Then there exists a spectral triple (A, 7i, D) 
with A = C" satisfying all the axioms, and such that the resulting distance on the set of pure 
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states of A is given by the numbers dij. 

To proceed, we shall first prove the following lemma. 

Lemma 14. There is a spectral triple (A, 7i, D) with A = C" satisfying all the axioms such 
that 

||[D,7r(a)]|| = sup (18) 
where a = (a\, . . . , a n ) G C n and n denotes the representation of A on TC. 

Proof. The proof is by induction on n. 

The first non trivial case is n = 2. We take A 2 = C 2 and 7i 2 = C 3 . The representation 
7r 2 and the chirality \2 are both diagonal and are given by 7r 2 (a;i,x 2 ) = diag(xi, x 2 , x 2 ) and 
X2 = diag(l, —1, 1). The Dirac operator D 2 and the charge conjugation J 2 are defined as 

V 2 = 

where C is the complex conjugation and we set = if d 12 = oo. 

In the finite case, all axioms reduce to the reality, first order, orientability and Poincare 
duality axioms^. In the present case, the first two are commutation relations easy to check due 
to the commutative nature of the algebra. The orientability axiom is fullfiled by writing the 
chirality as 

X 2 =7T 2 (l,-l) l 7 2 7T 2 (-l,l) l 7 2 - 1 . 

The multiplicity matrix is 
which is non degenerate and thus Poincare duality holds. Finally, one easily checks that 















di2 


|. H 
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:) 
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V i 








1 

1 -1 



[D 2 ,7T 2 (X)} 



Fi - x 2 \ 



Let us now assume that {A n ,7i n , D n ) together with 7r n , \n and J n have been constructed 
for n > 2. To build the order n + 1 spectral triple, we merely imitate the n = 2 construction. 
Let us take A = C n+1 and 
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with Tt l n = C 3 , Vz, n. With respect to the previous decomposition all operators O corresponding 
to V, 7T, x an d J are block diagonal and defined inductively as 

a = a-i®(Vo; 

As in the n = 2 case, we define 

7T l n (xi, x n ) = diag(xj, x n , x n ) , Xn = diag(l, -1, 1). 
The Dirac operator D n and the charge conjugation J n are defined as 















I) 

din 1 


V = 
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/ 












Then it is easy to check that all axioms but Poincare duality hold using the induction assumption 
and the block diagonal nature of the construction. 
The multiplicity matrix of this spectral triple is 

/0 1 ... 1 \ 

i -l '•. ; 
; ■-. ■-. i 

\ 1 ... 1 -(n-l) j 
If iV is any positive integer, we can always consider the trivial spectral triple (C n , C n © C^, 0) 
with obvious representation and charge conjugation and whose chirality is equal to -1. If we 
take the direct sum of this spectral triple with (A n , 7~Cn, D n ), the resulting multiplicity matrix 
is fi n + NI n , which is non degenerate for iV sufficiently large. Accordingly, Poincare duality 
will be satisfied. 

Finally, the computation of the norm of the commutator || [D n , TT n (a)] || follows easily from 
the block diagonal structure and the induction assumption. ■ 

To complete the proof of proposition 13, we use the previous lemma to construct a spectral 
triple (A,7i,D) fulfilling condition (Tl8|). If a verifies the norm condition, then — aj\ < dij, 



so d(i,j) < d^. 

Furthermore, if we fix any two points such that d^ < oo and take Xk = &%}. (which is finite 
thanks to triangular inequality) one has |xj — Xj \ = dy. An other application of the triangular 
inequality yields \dik — du\ < d^i for any k and I so that || [D, 7r(x)] || < 1 by (|T8D. This shows 
that < d(i,j) so that the equality holds when dy is finite. 

If d^ is infinite, so are dik and djk for any k. Thus, the inequality fllTf ) does not constraint 
Xi and Xj since the corresponding matrix element of D vanish. Therefore, we can send |xj — Xj\ 
to infinity and we also have d(i,j) = d^. ■ 
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X Conclusion. 



As a conclusion of previous discussion, we may say that once given A = C™, there is no 
constraint arising from the axioms of noncommutative geometry. Such constraints may only 
appear if one imposes some extra conditions, such as fixing 7i = C n as we did in the discussion 
of the three and four point cases. We stress that we only showed that the map which associates 
a metric to a Dirac operator is surjective. In a discrete analogue of a quantum theory of gravity 
based on eigenvalues of the Dirac operatorsEl, one also needs to know how many Dirac operators 
correspond to a given metric, as well as the possible relations between their spectra. 

A naive question remains unanswered: what does these distances mean? According to 
spectral action principle^, Dirac operator encodes both physics and metrics. In the standard 
model, computation of the action leads to the lagrangien of the full standard model and the 
Einstein-Hilbert action of general relativity. So the coding of physics makes sense. But what 
about the coding of metrics ? For the time beeing, the answer is clear in the simple continuous 
case where noncommutative distance is the geodesic distance. Distance in the geometry of 
standard model should have physical meaning, but it has not been explicitly computed yet. 
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Appendix. 



Coefficients of V e ff(x,y) in the general four point case. 



V 4 (x) 
V 3 (x) 



4(d 3 d 4 - d 2 d 5 ) 2 (d 4 2 dl + d 2 2 (d 4 2 + dp) 



d"2 ^3 d ^ d^ d§ 



8x(d 2 d 5 — d 3 d i )(d ? ,d A d 5 d & (d\ + df) + di(d 2 d 3 d 4 (dl + dp — d\d$d\ — d\d^{d\ + 2d§))) 

^1 ^2 ^3 d^d^d^ 
Ax 2 

«1«2«3«4«5«6 

-2d x d A d$(d 2 d A d?>(4 - 2^) + d 3 <^<^ + 44(^4 + 3d§)) 

0^4(0^3^4 - + 4(4(4 + 4 + 4) - Qd 2 d 3 d 4 d 5 + 0^2(^4 + Qdl)))} 

A{d 2 {d 2 {d 2 + d 2 + dp + d\{d\ + 2rf|)) - 2d 2 d 3 d A d 5 d 2 + + <g(<jjj + 2d 2 ))) 

d 2 d\d\d 2 d 2 

8x 3 (di4 - dsd^jd^dsdjjd 2 + dp - {d\d 2 d\ - d 3 d 4 (d 2 + d 2 )d 5 + d 2 {2d\ + d 2 )d 2 )d 6 ) 

did 2 d 2 d\ d^djQ 

MWsWj - ^5) + d^dljd 2 + g) + rf§(rf 2 (rfi + 14) ~ rf 3<W)) 

C?l C?2 ^3 ^4 ^6 

T / / \ _ . / .-2 , J-2 , . -2N , 4x 4 (rf^ + g(g + d 2 )){d 3 d 4 - rfx4) 2 

Vo(s) - 4(rf 3 +rf 6 +d % ) + rfi^p2 

Ax 2 (d\(2d 2 d 2 + c^l + rig)) - 2d 1 d 3 d A d 2 d 6 + g(g(g + 2d 2 ) + g(g + j§ + 4))) 

d~\_ d^ d^ d^ d^ 

Computation of d(l, 2) when ^ = ^ = 00. 



V^x) 



+ 




, , , , . , / '•':; "!: , i<7;l7;i " d\'lf,) \d\ + d\ 



4 + ^i ' V(d 3 rf4-did6) 2 V d 3 + ^' 

di|d 3 - di\ ± |d 4 (di + d e )\ d 3 (di + d % ) 

TJ2 = — ; =■ , z 2 = ±- 



\J(d 3 -d A f + {d 1 + d & ) 2 ' sJ(d 3 -d A f + {d 1 + d & f 
_ di\d 3 + di\ ± \di(di - d 6 )\ _ d 3 (d 1 - d % ) 

V3 — ; = , Z 3 



^J(d 3 + d 4 ) 2 + (d 1 -d 6 ) 2 y / (rf 3 + 4) 2 + (rfi-4) 2 

The choice of the signs depends on the sign of expressions in modules. 
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Computation of d(l, 3) when j- = j- = oo. 



d\d§\/ d\ + dl d\ 
— — ; — ; ; — ; — j z o 



Z2 



z 3 



did e - d 3 d 4 " yjdl + dl ' 
d\ d 3 d 4 ^J d\ -\- d\ 

\fd\ + d\ ' d 3 d 4 - did 6 

■ i a a , a a \ d l (d 3 + d A ) 

x 2 = sign(dia 3 + d 4 d 6 ) 



\J (d 3 + d 4 ) 2 + (di - d 6 ) 2 

d 3 (d 3 ^ (did 3 + d 4 d 6 ) 2 ± d 6 (d 3 (d 3 + d 4 ) - d x d§ + d 6 2 )) 

(d 3 + d 4 f + (di - d 6 ) 2 (d 3 2 + d 6 2 ) 
di(d 3 — d 4 ) 



x 3 = sign(did 3 + d 4 d 6 ) 



\J{d 3 -d 4 ) 2 + {d 1 + d & ) 2 
d?,{d 4 \J (d~id 3 + d 4 d 6 ) 2 ± d 6 (d 4 (d 4 - d 3 ) + d 1 (d 1 + d 6 ))) 



(d 3 d 4 - did 6 )\J (d 3 - d 4 ) 2 + (d 1 + d e ) 2 
The choice of the signs depends on the sign of expressions in modules. 
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